Tuning Kinetic Magnetism of Strongly Correlated Electrons via Staggered Flux 



O 

o 



Yi-Fei Wang^'^ Chang-De Gong^^^ and Z. D. Wang^'i 
^National Laboratory of Solid State Microstructures and Department of Ptiysics, Nanjing University, Nanjing 210093, China 
^Department of Physics and Center of Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
(Dated: February 1, 2008) 

We explore the kinetic magnetism of the infinite- C/ repulsive Hubbard models at low hole densities 
on various lattices with nearest-neighbor hopping integrals modulated by a staggered magnetic flux 
Tuning (j) from to tt makes the ground state (GS) change from a Nagaoka-type ferromagnetic 
state to a Haerter-Shastry-type antiferromagnetic state at a critical <j!>c, with both states being of 
kinetic origin. Intra-plaquette spin correlation, as well as the GS energy, signals such a quantum 
criticality. This tunable kinetic magnetism is generic, and appears in chains, ladders and two- 
dimensional lattices with squares or triangles as elementary constituents. 
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Introduction. — The study of magnetism may begin 
with a Heisenberg-type model in which localized spins 
interact with each other. However, the fundamental in- 
teraction between electrons is the Coulomb-like repul- 
sion, while the Pauli principle could induce a spin de- 
pendence when electrons are not localized. The Nagaoka 
theorem provides us a rigorous mechanism that the sat- 
urated ferromagnetic (FM) state of kinetic origin is the 
unique ground state (GS) when a single hole is inserted 
into the half-filled Hubbard model with an infinite on-site 
repulsion J7 [1,] . In the Nagaoka's problem, the sign of the 
hopping amplitudes around the smallest closed loop in a 
lattice, (Sloop, is necessarily positive. 

In a recent letter, Haerter and Shastry (HS) have made 
important progress in an opposite situation to Nagaoka's 
problem, a single hole moving in the infinite-C/ Hubbard 
model on two-dimensional (2D) triangular lattices with 
frustrated hopping, i.e., iSioop of an elementary triangle 
is negative Q . Through numerical exact diagonalization 
(ED) studies of finite lattices, they found that the motion 
of a single hole with the electronic frustration leads to 
weak metallic antiferromagnetism (AF) of kinetic origin. 

There are still two important issues to be resolved. One 
is whether there is a tunable transition between the two 
opposite kinetic magnetisms. The other is whether or not 
the two kinetic magnetisms can appear in the systems 
with more holes or with finite hole densities, which are 
much closer to experimental realities. Therefore, we are 
motivated to study the infinite-C/ Hubbard models at low 
hole densities on various lattices with nearest-neighbor 
hopping integrals modulated by a staggered magnetic 
flux ±0. In these systems, 5ioop = exp(i0) can be tuned 
from +1 to —1 via the Aharonov-Bohm effect, (p = and 
(j> = TT correspond respectively to Nagaoka's and HS's 
problems, and the spatial periodicity is preserved when 
(j) varies from to tt. 

Our motivations also come from three other aspects: 
our problem is relevant to future experiments applying 
an external magnetic field on artificial lattices of quan- 




FIG. 1: (color online), (left) Lattice structures with squares 
as elementary plaquettes: diamond chain, square ladder and 
2D square lattice, (right) Lattice structures with triangles as 
elementary plaquettes: sawtooth chain, trestle ladder and 2D 
triangular lattice. A plaquette filled with color has a flux (f), 
and a blank one has a flux — (ji. Each arrow represents a phase 
shift in its direction. 



turn dots Q or creating an artificial magnetic field on 
optical lattices of ultracold fermionic atoms an effec- 
tive magnetic flux can be induced intrinsically through 
some mechanisms in strongly correlated systems @, Q; 
and our problem presents a concrete example to study 
quantum phase transitions in both finite and infinite sys- 
tems with a tunable parameter 0. 

On the basis of ED calculations of finite systems and 
analytical estimations of infinite systems, we find that 
tuning (f) from to tt makes the GS change from a 
Nagaoka-type FM state to an HS-type AF state at a criti- 
cal 0c- Intra-plaquette spin correlation, as well as the GS 
energy, signals such a quantum criticality. This tunable 
kinetic magnetism is generic, and appears in many lattice 
structures with squares or triangles as elementary con- 
stituents: elementary square and triangle, diamond and 
sawtooth chains, square and trestle ladders, 2D square 
and triangular lattices. 
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TABLE I: Elementary square and triangle: GS properties of 
the cases with Nh ~ 1. 

Elementary Square 



chain, E^\^{(l)) 
chain, E^M{(f>) - 



— 2\/2 cos((/)/4); while for a sawtooth 
-(l + \/5) cos{(t>/3). 





Stot of GS 


{Si ■ Sj)n,n. 


GS Energy 


<<?!)< 7r/3 


3/2 


1/8 


-2cos((^/4) 


7r/3 < (f> < TV 


1/2 


-1/8 


-2cos(0/4 - 7r/6) 


Elementary Triangle 


<P 


Stot of GS 


{Si ■ Sj)n,n. 


GS Energy 


0<<j}<n/2 


1 


1/12 


-2cos((^/3) 


tt/2 < 4> <tv 





-1/4 


-2cos(</>/3 - 7r/3) 



Model Hamiltonian. — The infinite-U Hubbard model 
with a staggered flux (SF) can be written as: 



H 



t ^ e''^-(l-n,._,)c],c,,(l-n,,_,)+H.c. (1) 



where the hopping integral t is positive and is taken as 
the unit of energy, Cio- (c|^) is an electron annihilation 
(creation) operator on site i of spin cr =t or J,, is the 
electron number operator, and {ij) refers to two nearest 
neighboring sites. The magnetic flux per plaquette (the 
summation of along four links around a plaquette) 
is given by ±0 alternatively in neighboring plaquettes, 
with (j) in units of (/)o/27r (</)o = hc/e is the flux quantum). 
Since the system is symmetric under the transformations 
(j) —f —(f) and — > 27r — it is sufRcient to restrict (p in 
the interval [0,7r]. The hole number is = Nl — N^, 
where and Ng are the numbers of sites and electrons, 
respectively; the hole density is denoted by a:; = Nh/NL. 

The symmetric gauge is chosen for this SF as shown 
in Fig. [U and the corresponding periodical boundary 
conditions are adopted. It should be noted that trans- 
lational and macroscopic time-reversal symmetries make 
our model distinct from the uniform flux case, where the 
energy spectrum exhibits the fractal Hofstadter butterfly 
and the uniform flux induces a saturated FM from sta- 
tistical transmutation at a flux quantum per electron Q . 

Toy models of elementary square and triangle. — As an 
illustration of basic physics, we first consider toy mod- 
els of elementary square and triangle with a single hole 
{Nh = 1) which can be solved analytically. For the 
elementary square, when (j) increases from to tt, at 
(j)c — 7r/3, the GS transits from a state with the maxi- 
mum total spin Stot = | (a Nagaoka FM) to a state with 
'S'tot = |, and the nearest-neighbor (n.n.) spin correla- 
tion changes from | to — |. For the elementary triangle, 
when (j) increases from to tt, at (j)c — 7r/2, the GS tran- 
sits from a ^tot = 1 state to a Stot — state, and the 
n.n. spin correlation changes from to —j. 

Diamond and sawtooth chains. — On a periodic lattice 
with an SF, the lowest kinetic energy of a single hole in a 
saturated FM spin background, £'p^(0), can be obtained 
via the Fourier transformation. For an infinite diamond 




FIG. 2: (color online). (a),(b) Diamond chains. (c),(d) Saw- 
tooth chains, m's and rescaled GS energies versus in the 
cases with various sizes and hole numbers. 

We will present typical ED data obtained by using the 
Spinpack package 9] and exploiting the translational in- 
variancc in the subspacc with fixed S^Qt (0 fo^' even N^, 
and i for odd N^). The commonly used quantity in eval- 
uating Nagaoka FM, i.e. the GS Stot, is strongly de- 
pendent upon the boundary conditions chosen and the 
even/odd parity of Nh [1, llO) Hi 12\- We therefore con- 
centrate on two derived quantities which are not sensitive 
to the boundary condition or the parity of N^. 

The first quantity, m, is used to measure the spin 
correlations intra a plaquette. For a square plaquette, 
mi = i + (5*2 • 53) -I- (5*2 • 54) + (S'a • 54) (with four clockwise 
sites 1, 2, 3 and 4, and (...) means the GS average) [3], 
and the average over four sites gives m\j = \ "^j- 
On the square lattice, mu = 1 in a classical FM state, 
mn = in a classical Neel AF state, while mu = \ 
if there is no spin correlation. For a triangular plaque- 
tte, TOi = ^ + 2(52 • S3), and the average over three 
sites gives — triangular lattice, 

toa = 1 in a classical FM state, = | in a classical 3- 
sublattice 120° AF state, while = 5 if there is no spin 
correlation. The other quantity, the rescaled GS energy 
E{(j))/\E{0)\, is used to compare the non-analyticities in 
GS energies of various cases. 

For the 12-site and I8-site diamond chains with re- 
spectively Nfi = 1 — 2 and Nh — 1 — 4, when 
changes from to tt, toq's drop almost abruptly near 
a (/)c « 7r/3 [Fig. H^a)]. (Note that 0c = n/3 for the 
elementary square with Nh — 1.) Meanwhile, the GS 
energies also show clearly certain non-analyticities near 
tt/S [Fig.EKb)]. 

For the 12-site and 18-site sawtooth chains with re- 
spectively Nh = 1 — 3 and Nh = 1 — 4, as seen from 
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Figs. [UJc) and (d), the abrupt drops of t^a's and non- 
analyticities in GS energies occur near a, (j>c ~ 7r/2. (Note 
that (j)c — 7i'/2 for the elementary triangle with Nf^ = 1.) 

In the four cases of diamond chains {Nl = 12 with 
Nh = 1, Nl = 18 with Nh = 1-3) and all seven 
cases of sawtooth chains, m's change very little (and the 
systems possess intra-plaquette FM correlations since all 
m's satisfy m > ^) when cj) varies from to (f)c, and 
the curves of rescaled GS energies are very close to that 
of E^^^{(j3)/\E^\j^{0)\ [the continuous curves in Figs.HJb) 
and (d)]. (The ED data also tell us that in each case 
with odd Nfi, the GS Stot takes the maximum value 
iVe/2 when < ^ (pc-) Since (j> ~ corresponds 
to the Nagaoka FM (5ioop — +1) and m's change lit- 
tle for < (/) < 0c, we extend the notion of Nagaoka 
FM here and such a state (with < < (j>c) is referred 
to as a Nagaoka-type FM. When (p varies from (pc to 
TT, all to's show obvious intra-plaquette AF correlations 
{rriu < J for diamond chains, and < 5 for sawtooth 
chains). Since cp — tt corresponds to the HS's problem 
(^loop = ^1) and m's change little as (/)c < < tt, we 
also extend the notion of HS AF and such a state (with 
(pc < (p < tt) is referred to as an HS-type AF. 

Square and trestle ladders. — For the single-hole Na- 
gaoka FM in an infinite square ladder with an SF, 
+ A cos(0/2); while in an infinite tres- 
(0) = -4cos(0/3). 
For the 8 x 2, 10 x 2 and 12x2 square ladders with re- 
spectively Nk. — Nh — 2 — 4 and Nh — l — 2,as seen 
from Figs, ^a,) and (b), the abrupt drops of mg's and 
non-analyticities in GS energies occur near different (/>c's, 
with (pc versus hole density x varying rather smoothly. 

For the 6x2 and 9x2 trestle ladders with respectively 
TV/, = 1-3 and A^?, = 1-4 [Figs.E^c) and (d)], the abrupt 
drops of mA's and non-analyticities in GS energies also 
occur near different (pc's, with still a smooth curve of (pc 
versus x. 

In these cases of square and trestle ladders, m's change 
very little and satisfy m > ^ when (p varies from to 0c , 
the curves of rescaled GS energy are very close to that of 
'^fm(^)/I-^fm(0)Ii and these states are of Nagaoka-type 
FMs; when (p varies from (pc to tt, all m's show obvious 
intra-plaquette AF correlations, and these states are of 
HS-type AFs. 

It is tempting to estimate (pc in the limit of low hole 
density {x — ^ 0). Such a task can be partly fulfilled 
with the Brinkman-Rice (BR) approximation lj| . For a 
single hole in an infinite Neel AF spin background, the 
BR approximation accounts the dominant contributions 
to the self energy of single- hole Green's function, i.e., 
the retraceable paths without any closed loop. Such an 
approximation leads to a hole band edge (i.e., the lowest 
single-hole kinetic energy) E^p = —2\Jz — 1, where z is 
the coordination number. 

Through comparison between E^ and i?p^((/>) of a 
single hole in infinite ladders, one can obtain a rough 




FIG. 3: (color online). (a),(b) Square ladders. (c),(d) Trestle 
ladders, m's and rescaled GS energies versus (p in the cases 
with various sizes and hole numbers. The insets in (b) and 
(d) show 0c vs. hole density x. 



0. For an infinite 



^fm('/') gives a 



estimation of (pc in the limit x 
square ladder, z = 3, and E\ 

0jf^ = 2arccos(3/4) w 0.467r. While for an infinite tres- 
tle ladder, z = 4, and E^ = £^fm(0) gives a 0|?^ = 7r/2. 

Square and triangular lattices. — For the single- hole 
Nagaoka FM in an infinite square lattice with an SF, 
£'p^(0) = — 4cos(0/4); while for an infinite triangular 
lattice, £^fm(<^) = -6cos(0/3). 

For the 4x4 and 6x4 square lattices with respectively 
Nh = 1-4 and TV,, = 1-2 [Figs.^a) and (b)], or the 3x3 
and 6x3 triangular lattices with respectively Nh — 1 and 
Nh~l — i [Figs.HKc) and (d)], the abrupt drops of m's 
and non-analyticities (or changes of concavities) in GS 
energies also occur near different 0c 's. In all these cases, 
the GSs are of Nagaoka-type FMs for < < 0c, and 
are of HS-type AFs for 0c ^ < tt. The curve of 0c vs. 
X of the square lattices [the inset in Fig. [4Kb)] approaches 
a critical doping Xc ~ 0.3 at 0c = 0, which agrees well 
with the quantum Monte Carlo studies on the instability 
of saturated Nagaoka FM against doping [12]. 

A rough estimation of 0c in the limit x ^ with the 
aid of the BR approximation is: 0^^ — 27r/3 for an in- 
finite square lattice, and 0^^ = 3 arccos(-\/5/3) ~ O.TOtt 
for an infinite triangular lattice. 

Long-range spin correlations and ordering. — As seen 
from the above, the two quantities, m and £'(0)/l£'(O)l, 
describe well the transitions from the Nagaoka-type FM 
to the HS-type AF. The ED data of the six kinds of 
lattices also tell us that for each case with odd Nh, the GS 
S'tot always takes the maximum value Nc/2 and {Si-Si+r) 
is positive and almost a constant for any range r when 
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FIG. 4: (color online). (a),(b) Square lattices. (c),(d) Trian- 
gular lattices, m's and rescaled GS energies versus <j) in the 
cases with various sizes and hole numbers. The insets of (b) 
and (d) show <j!>c vs. hole density x. 



< <C 0c, namely, the Naogaka-type states are long- 
range ordered FM. Now we would take a closer look at 
the long-range spin correlations of the cases with even 
NhS. We focus on the three cases of square ladders (in 
which there are the longest-range spin correlations) as 
examples: A^l = 10 x 2 and Nh = 2, 4, A^i, = 12 x 2 and 
Nh = 2. 



but extended (in the whole lattice) AF correlations. 

Summary and discussion. — As for the infinite-C7 repul- 
sive Hubbard models on the considered various lattices at 
low hole densities {x < 0.3) with n.n. hopping integrals 
modulated by a SF ±ip, ED calculations of finite systems 
and analytical estimations of infinite ones enable us to 
find: (i) tuning from to tt makes the GS change from 
a Nagaoka-type FM to an HS-type AF at a critical (pc, 
and both states are of kinetic origin; (ii) this tunable ki- 
netic magnetism is generic, regardless of the parity of N^, 
and appears in various lattice structures with squares or 
triangles as elementary constituents, i.e. elementary pla- 
quettes, chains, ladders and 2D lattices; (iii) near a 4>c, 
the intra-plaquette spin correlation drops abruptly from 
FM to AF, and the GS energy exhibits non-analyticity; 
(iv) while the (jjcS of chains are insensitive to x, the (jjc's 
of ladders or lattices versus x exhibit a decreasing and 
almost smooth function; (v) the Nagaoka-type FM has 
long-range FM ordering, while the HS-type AF has de- 
cayed but extended AF correlations. 

In an artificial lattice of quantum dots, a large lattice 
constant would enable us to observe this effect at a mod- 
est flux strength of a few Tesla Q. This effect could 
also be realized in optical lattices of ultracold atoms if 
appropriate phase factors are introduced for hopping in- 
tegrals by laser assisted hopping, lattice tilting and other 
experimentally accessible techniques 4] . 

This work was supported by the NSFC and RGC 
grants of Hong Kong. 




FIG. 5: (color online). Long-range intra-chain (x) and inter- 
chain (o) GS spin correlations of square ladders at <j!> = and 
(j) = TT. The insets in (b) and (f) enlarge the middle parts. 

At (j) = [Figs.[5i;a),(c),and (e)], the GSs of the three 
cases are singlets (^tot = 0), however, there is an evi- 
dence of long-range FM in {Si ■ Si+r) and each singlet 
GS actually consists of two FM domains with opposite 
magnetization as pointed out by Troyer et al. [ll[ . While 
at = TT [Figs. [njb),(d), and (f)], the GSs have decayed 
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